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We perform Raman scattering experiments on natural graphite in magnetic fields up to 45 T, 
observing a series of peaks due to interband electronic excitations over a much broader magnetic 
field range than previously reported. We also explore electron-phonon coupling in graphite via 
magnetophonon resonances. The Raman G peak shifts and splits as a function of magnetic field, 
due to the magnetically tuned coupling of the E2g optical phonons with the K- and if-point inter- 
Landau- level excitations. The analysis of the observed anticrossing behavior allows us to determine 
the electron-phonon coupling for both K- and iif-point carriers. In the highest field range (>35 T) 
the G peak narrows due to suppression of electron-phonon interaction. 

PACS numbers: 78.30.-j, 71.70.Di, 73.61. Cw, 76.40.+b, 78.20.Bh 



Electron-phonon coupling in graphene and graphite 
has been investigated for several years^^. The zone- 
centre, doubly degenerate, E2g phonon strongly interacts 
with electrons, resulting in renormalization of phonon 
frequencies and line broadening^^^. These are tun- 
able by electric and magnetic fields, through Fermi- 
energy shifts and Landau quantization. The Raman G 
peak is predicted to exhibit anticrossings when the E2g 
phonon energy matches the separation of two Landau 
levels (LLs). Both intraband (i.e., cyclotron resonance- 
like) and interband (i.e., magnetoexcitonic) transitions 
are allowed both in single-layer graphene (SLG)^*^ and 
bilayer graphene (BLG)^^. Interband magnetophonon 
resonance (MPR) has indeed been observed in magneto- 
Raman scattering on SLG on the surface of graphit o^^i^^ 
and non-Bernal stacked multilayer graphene on SiG^^. 

Graphite, a semimetal containing both electrons and 
holes even at zero temperature, is expected to exhibit 
even richer carrier-phonon coupling phenomena. Indeed, 
Ref. [13 recently reported magneto- Raman measurements 
on graphite up to 28 T, and observed inter-LL tran- 
sitions and signatures of MPR. As described via the 
Slonczewski-Weiss-McGlure (SWM) modeU^i — , graphite 
has a linear ("massless") dispersion for the hole pocket 
around the H point of the Brillouin Zone and a parabolic 
( "massive" ) dispersion for the electron pocket around the 
K point. Angle-resolved photoemission measurements 
provided evidence of such massless and massive quasi- 
particles in graphite^i. Near these high symmetry points, 
graphite's band structure can be approximated as a com- 
bination of SLG, describing the i^-point massless holes, 
and BLG, describing the i^-point massive electrons^^. 

Here, we report low-temperature magneto- Raman 
measurements of natural graphite in a magnetic field {B) 
up to 45 T, a range of fields much broader than any pre- 
vious study, to the best of our knowledge. We demon- 



strate a rich picture of MPR effects caused by coupling 
of the E2g phonon to both i^-point (SLG-like) and K- 
point (BLG-like) interband excitations. We also observe 
a series of electronic Raman excitations (i.e., emission of 
electron- hole pairs instead of phonons), including transi- 
tions involving the lowest, electron-hole mixed, LLs. We 
explain the entire, complex set of Raman-active inter- 
band excitations within a SWM approach. Furthermore, 
through quantitative analysis of the observed anticross- 
ing behaviors, we determine the strengths of electron- 
phonon coupling (EPG) for both i^-point holes and K- 
point electrons. Finally, in the highest magnetic-field 
range (>35 T), where all transition energies are far away 
from the E2g phonon energy, the G peak narrows, due to 
suppression of the EPC contribution to the linewidth. 

Raman spectra were collected on natural graphite 
(NGS Naturgraphit GmbH) in a backscattering geom- 
etry, with B up to 45 T [see Fig. 1(a)]. A 532-nm laser 
is coupled via an optical fiber to the low-temperature 
probe, and focused to a spot of < 20 /im, with a power 
of ~13 mW. The probe is inserted into a helium cryostat 
and placed in a 31-T resistive magnet or 45-T hybrid 
magnet. Under laser illumination, the temperature of 
the sample is stabilized at ~10 K. The unpolarized Stokes 
component of the scattered light is directed into the col- 
lection fiber and guided to a spectrometer equipped with 
a charge-coupled-device camera. Most of the data were 
collected with a spectral resolution of ~3.4 cm~^. How- 
ever, we used a spectral resolution r^{).b cm~^ to accu- 
rately measure the full width at half maximum of the 
G peak, FWHM(G), at selected magnetic fields between 
32 T and 45 T. Raw data contains the signal of interest 
from the sample in a smooth background coming from 
the fibers. At frequencies >1300 cm~^, the background 
is featureless and much smaller than the signal from the 
sample. We performed numerous tests to characterize 
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FIG. 1. (Color online) (a) Experimental configuration, (b) 
Unpolarized, background-corrected Raman spectrum of nat- 
ural graphite at zero magnetic field and 10 K. (c) Magneto- 
Raman spectra at 10 K for various magnetic fields. A series 
of interband electronic transitions are observed, as indicated, 
as well as a large G peak~1580 cm~^. Gray arrows point to 
weak low- frequency scattering peaks. 



the field and temperature dependence of the background, 
in particular at frequencies <900 cm~^, where parasitic 
scattering in fibers becomes comparable or higher than 
the signal from the sample. As the magnetic field in- 
creases up to 45T, we observe small broad-band changes 
reaching ±1% at ~650 cm~^. However, these could be 
due to other factors, such as long-term variation of the 
laser power or temperature drifts. Thus, we assume the 
background to be field independent, at least within our 
signal-to-noise ratio and use zero-field reference spec- 
tra to remove spurious signals due to scattering in the 
fibers. Decoupled SLG may exist on the surface of bulk 
graphitosis. To avoid contributions from such SLGs, 
we recorded the spectra of several locations to select a 
region with a bulk graphitelike Raman spectrum, as in- 
dicated by the 2D Raman peak shape^j^. This approach 
is opposite to that of Refs. il^ andllsL where the samples 
were scanned to find SLG Raman signatures. 

At 5 = 0, the first-order Raman spectrum of graphite 
is dominated by the G peak at ~1580 cm~^, as shown 
in Fig. 1(b), due to scattering by the doubly degenerate 
zone-center E2g phono n^S. As B increases, a number 
of peaks emerge, as shown in Fig. 1(c). These peaks be- 
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FIG. 2. (Color online) Intensity map of magnetic- field- 
dependent Raman intensity. Dotted lines represent the cal- 
culated energies of Raman- active, hn en, electronic exci- 
tations near the K point. Dashed lines represent asymmet- 
ric, hn e{n — 1) and hn e{n -\- 1) K point excitations 
coupled to the E2g phonon. Open circles indicate the weak 
low- frequency peaks labeled by gray arrows in Fig. 1(c). 



come sharper and move towards higher frequencies with 
increasing B. Similar peaks were previously reported in 
Ref. where Raman scattering of bulk graphite was 
measured in magnetic fields up to 6.5 T, but assigned to 
LLs in BLG. 

Figure 2 displays a set of spectra taken at 10 K as 
a function of 5 up to 45 T. The observed nearly linear 
B dependence suggests these features to be related to 
inter-LL excitations of massive carriers in the vicinity of 
the K point. The most intense peaks are attributed to 
the so-called "symmetric" inter-LL excitations, hn -^en 
or {n^n)^ i.e., the transitions from the nth hole to the 
nth electron LL^. Indeed, Ref. HH and HqI theoretically 
showed that symmetric inter-LL excitations are Raman 
active in both SLG and BLG. These symmetric transi- 
tions were previously observed and analyzed through an 
effective BLG model^^. 

In addition, we detect two extra electronic features 
below the (1,1) transition, indicated by open circles in 
Fig. 2. They are resolved at 45 T, as shown by gray ar- 
rows in Fig. 1, although their intensity is less than 10% of 
the (1,1) peak. We attribute them to the lowest inter-LL 
transitions, (1,0) and (—1,1), at the point. They can 
be considered as a special case of the weak lowest-energy 
Raman-active transition in BLG predicted in Ref. [29|. 

To validate our peak assignments, we calculate the en- 
ergies of interband, inter-LL transitions within the SWM 
model. This has seven tight-binding parameters, 70 to 75 
and A. Despite its extensive use over the past 50 years, 
the precise values of these parameters are still under de- 
bate. Without the trigonal warping effect represented by 
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TABLE I. SWM band parameters (in eV) extracted from re- 
sults in Fig. 2, in comparison with previously reported values. 

This work Ref . 30 Ref . 31 Ref . 17 Ref . 32 

70 3.06 (1) 3.1 3.18 (3) 3.08 (1) 3.16 (5) 

71 0.370 (5) 0.39 0.38 (1) 0.380 (2) 0.39 (1) 

72 -0.028 (4) -0.028 (4) -0.02 -0.020 (2) 

73 0.33 (1) 0.315 - 0.315 (1) 0.315 (15) 

74 0.080 (5) 0.041 (10) 0.08 (3) 0.044 (5) 0.044 (24) 
A + 275 0.130 (3) 0.15 (3) 0.064 (3) 0.084 (7) 

^ A + 275 - 272 = 0.22 (1) 



73, each LL can be obtained through a 4 x 4 Hamiltonian, 
which can be diagonalized for each n. Adding the 73 term 
mixes different LLs with indices n and n ± 3, making the 
dimension of the Hamiltonian infinite. We numerically 
calculate the LL energies by truncating this Hamiltonian 
into a finite ^ 400 x 400 matrix. We note that, for ma- 
trix sizes larger than 100 x 100, the gaps between energy 
levels at B =10 T change less than 0.1 cm~^. For higher 
magnetic fields, the results converges even faster. We ob- 
tained 70 from the position of the i7-point MPR and used 
the SWM parameters from Ref. |30| as the initial guesses 
for our fitting. To reduce the number of parameters, we 
fixed 70 and 72 and varied the others to fit the data. 

Table I compares our results with values ex- 
tracted from magnetotransport experiments^, in- 
frared magnetoreflectance spectroscopy^, magneto- 
Raman measurements^^, as well as values deduced from 
earlier infrared magneto-spectroscopy experiment^. 
Though the tight-binding parameters are not signifi- 
cantly different, our spectroscopic observation of both 
symmetric and asymmetric transitions, including the 
low-energy transitions involving the electron-hole mixed 
— 1 and LLs, enables an accurate determination of the 
SWM parameters. 

Close examination of the G peak in Fig. 2 reveals peak 
position modulations as a function of B. At a certain 5, 
the resonance condition En^n' = is met, where £^n,n' 
is the (n,n') transition energy and ftr is the E2g phonon 
frequency, and the phonon is "dressed" by the electronic 
transitioi>ii — . This coupling manifests itself as a series 
of avoided crossing^—. Specifically, the E2g phonon 
is allowed to couple with an (n^n') transition only when 
|n| — |n'| = ±1. To examine the data more closely, we 
fit the G peak with Lorentzians and plot the extracted 
peak positions and the second derivative Raman inten- 
sity in Fig. 3(a). The data reveals anticrossings at 34, 
31, 21, and 19 T, corresponding to the (2,1), (1,2), (3,2), 
and (2,3) transitions, respectively. At lower fields, the 
doublet structure due to the (3,4) and (4,3) transitions is 
smeared out and appears as a weak modulation of the G 
peak. Note that, when the symmetric (n,n) peaks cross 
the G peak, they appear unchanged, indicating the ab- 
sence of coupling. Furthermore, the central position of 
the G peak is also B dependent, exhibiting a modulation 
and broadening at ~30 T (Fig. 4), which we interpret as a 




FIG. 3. Magnetophonon resonance at the ii^ point, (a) Second 
derivative Raman intensity map. Circles: Positions extracted 
from Lorentzian fits. Dashed lines: Calculated energies of 
coupled i^^-electron-phonon modes, (b) G peak measured at 
different with multiple-Lorentzian fits. Inset: schematic 
energy diagram and MPR-coupled inter-LL transitions at the 
K point. 



signature of MPR of the asymmetric /il ^ and hi ^ —1 
iJ-point excitation with the E2g phonon. Finally, above 
35 T, where the decay of E2g phonons into electron-hole 
pairs is quenched by Landau quantization and electron- 
phonon interaction is suppressed, the G peak narrows 
to -4.4 cm-^ Our high-field value FWHM(G) is about 
twice the phonon- lifetime-limited linewidth at 5 = 0, 
^ph-ph _ 2.5 cm-i (Refs.lHandlai), indicating the pres- 
ence of another, probably disorder-induced, broadening 
mechanism. 

To analyze the observed MPR, we first focus on the 
G-peak sidebands, corresponding to coupled electron- 
phonon modes associated with iT-point electron asym- 
metric transitions [Fig. 3(a)]. The doublet anticross- 
ings at 34 and 31 T, corresponding to the (2,1) and 
(1,2) transitions, respectively, is most accurately resolved 
[Fig. 3(b)], and therefore, niost suitable for quantitative 
analysis. Following Refs. [12I and [13, we analyze the data 
via a two-coupled-mode model. 



^± = ^^^^^±\/(^^4^)' + i;^ (1) 

where Eq hftr - ^7r/2, £^n,n' ^^n,n' - ijn,n'/^, 7r 
(7n,nO is FWHM(G) [(n,n') transition], and g is the cou- 
pling parameter. Expressing the magnetic energy Jiwb at 
the K point within an effective BLG model, the coupling 
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FIG. 4. Magnetophonon resonance at the H point, (a) Peak 
position and (b) FWHM of the C-peak central Lorentzian 
component as a function of magnetic field. The model re- 
sults of MPR-induced modulation are shown by solid lines. 
Inset: schematic energy diagram and MPR-coupled inter-LL 
transitions at the H point. 

parameter g is given by: 

where a (= 2.46 A) is the graphite lattice constant, 70 
and 71 are tight binding parameters (see Table I), and 
Ib = ^JhJeB is the magnetic length. The dimensionless 
EPC Ar is defined following the notation of Refs. [lol and 

[all: 

^ Mh^r vl ^Mh^r 7^ 

where Au.c. is the graphene unit-cell area, M is the car- 
bon atomic mass, and vp = ^clIq = 0-99 x 10^ 
the Fermi velocity. {Dy^) is the deformation potential of 
the E2g phonon, which describes the modulation of the 
coupling energy 70 by C-C bond length variation. 

The position and linewidth of unperturbed phonons 
can be derived from the high-field (>35 T) spectra, i.e. 
Pos(G) = 1582.6 cm-i and FWHM(G) = 4.4 cm-^ We 
calculate the energies of asymmetric inter-LL transitions 
using the SWM parameters described previously. Fit- 
ting Eq. (1) to the anticrossings at 31 and 34 T yields 
7p = 44±6 cm-^ and g^^^ = 0.72±0.03 cm" VT- We can 
thus extract Ap^^ ^ 3.3 x 10~^, in excellent agreement 



with that previously derived from density functional the- 
ory, the zero-field FWHM(G), the doping dependence of 
Pos(G), and the slope of the phonon dispersions around 
F: X^!"^ ^ 3 X 10-U^^^ 

Finally, we analyze the 5-induced modulation of the 
central component of the G peak, shown in Fig. 4. The 
total peak-position modulation is ^6 cm~^, while the 
FWHM increases more than twice at ~30 T. This is con- 
sistent with MPR due to i^-point inter-LL transitions, 
(1,0) or (1,-1), assuming that the LL widths are larger 
than the coupling strength. The G-peak modulation at 
~20 T is a signature of the MPR effect involving (2,3) 
iC-point excitations. To deduce the EPC strength for the 
i^-point, we model the 30 T resonance with Eq. (1) using 
a SLG-like expression for g^^"^: 



The right-hand side of the resonance {B > 30 T) 
fits well with the model with 7r = 100 ± 10 cm~^, 
g^Q^^ = 3.2 ± 0.2 cm-VTi/2, and x\!^^ ^ 1.6 x 10"^. 
The discrepancy at lower fields is likely be due to the 
E2g renormalization via interaction with multiple inter- 
LL excitations, which cannot be spectrally resolved for 
B < 30 T. We note that Ap^^ is almost 20 times smaller 
than Ap^\ 

In summary, we performed high-field magneto-Raman 
experiments on graphite, observing strong magne- 
tophonon resonances. The G peak shifts and splits as 
a function of magnetic field as it sequentially resonates 
with certain electronic transitions. Analysis of the ob- 
served magnetophonon resonance effects allowed us to 
determine the strengths of electron-phonon coupling for 
both H- and i^-point carriers. The Slonzcewski-Weiss- 
McClure model provides an accurate description of all 
observed interband electronic excitations. In the highest 
field range (>35 T), the G peak narrows through reduced 
electron-phonon interaction. 
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